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I review recent studies that predict a realization of quantum liquid-crystalline orders in resonant 
atomic gases. As examples of such putative systems 1 will discuss an s-wave resonant imbalanced 
Fermi gas and a p-wave resonant Bose gas. In the former, the liquid-crystalline smectic, nematic and 
rich variety of other descendant states emerge from strongly quantum- and thermally- fluctuating 
Fulde-Ferrell and Larkin-Ovchinnikov states, driven by a competition between resonant pairing and 
Fermi-surface mismatch. In the latter, at intermediate detuning the p-wave resonant interaction 
generically drives Bose-condensation at a finite momentum, set by a competition between atomic 
kinetic energy and atom-molecule hybridization. Because of the underlying rotationally-invariant 
environment of the atomic gas trapped isotropically, the putative striped superfluid is a realization of 
a quantum superfluid smectic, that can melt into a variety of interesting phases, such as a quantum 
nematic. I will discuss the corresponding rich phase diagrams and transitions, as well the low-energy 
properties of the phases and fractional topological defects generic to striped superfluids and their 
fluctuation-driven descendants. 



I. INTRODUCTION 
A. Resonant atomic gases 

Experimental progress in trapping, cooling and co- 
herently manipulating Feshbach-resonant atomic gases 
opened unprecedented opportunities to study degener- 
ate strongly interacting quantum many-body systems in 
a broad range of previously unexplored regimes [THS]. 
These include paired fermionic superfluids (SF) [6Hl5]). 
the associated Bardeen-Cooper-Schrieffer (BCS) to Bose- 
Einstein condensation (BEC) crossover |31 [T2H21], Bose- 
Fermi mixtures [25]. bosonic molecular superfluids p6H29] . 
and many other states and regimes [30] under both equi- 
librium and nonequilibrium conditions |23l f?Tl I32j . 

Because degenerate atomic gases are free of the un- 
derlying crystalline matrix of the sohd-state materials 
(though one can be imposed through a highly tunable 
optical lattice potential[T]), among this rich variety of 
states, they admit phases that spontaneously partially 
break continuous spatial symmetries and thereby exhibit 
concomitant strongly fluctuating Goldstone modes with 
corresponding rich phenomenology. Resonant atomic 
gases are thus uniquely suited for a realization of quan- 
tum liquid-crystalline states of matter, that have been 
somewhat of a holy-grail dating back to their studies in 
solid state materials such as the striped states in high-Tc 
superconductors, nickelates and other strongly correlated 
doped Mott insulators ^33 38J, heavy-fermion and organic 
superconductors [Ml SO], spiral states in helimagnets [4ll- 
133] , and a two-dimensional electron gas with a partially- 
filled Landau level [lilH51H7] . In isotropic traps the pu- 
tative quantum liquid-crystal order is expected to exhibit 
all the complexity of fluctuations and topological defects 
of conventional (mesogenic) liquid crystals [3H1 BS]) but 
with the added enrichment of the accompanying quan- 
tum (off-diagonal) order of a superfluid. Another, not 
insignificant virtue is that (in contrast to other e.g., solid- 
state or nuclear matter systems) these dilute gases are 



extremely well-characterized at the two-body level, and 
are therefore described by microscopic (as opposed to ef- 
fective) Hamiltonians with well-known couplings. 



B. Candidate systems 

Recent theoretical studies have predicted a number 
of such quantum liquid crystal realizations in degener- 
ate atomic systems, that in addition to internal sym- 
metries partially break spatial svmmetries |50l - f5B] . These 
are typically driven by strong resonant and competing 
interaction that frustrates a spatially homogeneous and 
isotropic superfluidity. Known examples include bosonic 
and paired fermionic superfluids, where spatial order is 
driven by (i) dipolar interaction 60J , (ii) pseudo-spin- 
orbit interaction [57] (realized through hyperfine states 
coupled by Raman transitions[SH]), (iii) p-wave resonant 
interaction 1531 159] . and (iv) a Fermi surface mismatch 
(realized through species number and/or mass imbal- 
ance), that leads to the Larkin-Ovchinnikov- Fulde-Ferrell 
finite- momentum pairing [5(7ll51| . as well as their strongly 
fiuctuating descendent states. In this brief review, I will 
focus on the last two realizations, and will discuss the as- 
sociated microscopic models that I believe can realize a 
quantum superfluid liquid-crystal order, their phase be- 
havior, fluctuations, topological defects, and a variety 
of experimental predictions and signatures. Much of the 
discussed low-energy phenomenology is shared more gen- 
erally by systems exhibiting quantum liquid crystal or- 
ders in isotropic trap . For a more complete account, 
I refer the reader to the original literature and the more 
extensive reviews [21 [5]. 

I will not discuss the finite momentum states that 
depend on the lattice for their realization and stabil- 
ity, such as the p-band and FFLO superfiuids in opti- 
cal lattices 61 -65 . These are fascinating states, but are 
less relevant from the liquid-crystal perspective of this 
review. 
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II. IMBALANCED RESONANT FERMI GASES 
A. Background 

The most widely explored candidate for a realization 
of quantum liquid crystal order is a species-imbalanced 
Feshbach-resonant Fermi gas, [11121 IS], though it is only 
very recently that it was formulated and explored in 
these liquid-crystal terms [661167] , These studies build on 
well-explored two-species Feshbach-resonant Fermi gas, 
that exhibits paired superfluidity, that can be tuned be- 
tween a weakly-attractive Fermi-surface-driven BCS and 
a strongly-attractive molecular BEC superfluids[3l |4]. 
While at r = a balanced gas exhibits no qualitative 
change of state, a quantitatively accurate description of 
this crossover, particularly around the strongly interact- 
ing and universal unitary regime (where in a vacuum a 
two-particle bound state first forms and the s-wave scat- 
tering length diverges) has presented a considerable chal- 
lenge with much recent progress. 

A species-number (and mass) imbalance in the two 
atomic hyperfine-states mixture offered a new extremely 
fruitful experimental knob [68H71j . The imbalance frus- 
trates pairing j72l477j . driving quantum phase transitions 
out of the paired superfluid to a variety of possible ground 
states and thermodynamic phases [SH [75H55] . This 
rekindled considerable theoretical activity in the con- 
text of species-imbalanced resonant Fermi gases [Ml4106| . 
The corresponding imbalance versus detuning BCS-BEC 
phase diagram, illustrated in Figjl] (and its extension to 
finite temperature) is now well-established[79 l 11001 TlOSj . 
showing qualitative agreement with experiments |68fm] . 
More recently, considerable progress has been made to- 
ward establishing quantitative details of this phase di- 
agram through analvtical[107'- TTTl . numerical [11 2) 1113] 
and experimental approaches ^701 17T1 1114j . 

The identification of the number species imbalance 
with the magnetization of an electronic system, and the 
chemical potential difference with an effective Zeeman 
energy, connects these atomic gases studies with a large 
body of research on solid state electronic superconduc- 
tors under a Zeeman field 50, iSH I115H1171 . as well as 
extensively studied realizations in nuclear and particle 
physics [TTSHUI]- The obvious advantage of the newly- 
realized atomic system is the aforementioned tunabil- 
ity, disorder-free "samples", and absence of the orbital 
part of the magnetic field, that always accompanies a 
solid-state charged superconductor in a magnetic field. 
In these neutral paired superfluids the orbital field ef- 
fects can be independently controlled by a rotation of 
the atomic cloud |122] . 

As illustrated in Figjl] among many interesting 
features, such as the gapless imbalanced superfluid 
(SFM) f7WST\ ITUT]. ubiquitous phase separation [H Ei 
ITOOl ITOT] . tricritical point[IOOl [Ml ^06\, etc., ob- 
served experimentally [68U7I] and studied extensively 
theoretically [5]: the interaction-imbalance BEC-BCS 
phase diagram is also predicted [7S1 EH IH HOOl HMl 



I105j to exhibit the enigmatic Fulde-Ferrell-Larkin- 
Ovchinnikov state (FFLO)|5Ql[51]. First predicted in the 
context of solid-state superconductors over 45 years ago, 
the FFLO states has so far eluded a definitive observa- 
tion, though some promising solid state [123, and quasi- Id 
atomic [T24] candidate systems have recently been real- 
ized. 

At its most generic level the FFLO state is a fermionic 
superfluid, paired at a finite center of mass momentum. 
It spontaneously "breaks" gauge and translational sym- 
metry, a periodically-paired superfiuid (superconductor) , 
akin to a supersolid |125tn^5] , and thus can appropriately 
be called a pair-density wave {P'DW)^M USD]- This 
state can be equivalently thought of as a periodically or- 
dered micro-phase separation between the normal and 
paired states, that naturally replaces the macro-phase 
separation [711 1121| ubiquitously found in the BCS-BEC 
detuning-imbalance phase diagram [79 l 11001 HOlj I103j . 

Microscopically, it is driven by Fermi surface 
mismatch[SDl[5T] due to an imposed pairing species num- 
ber (and/or mass [132] ') imbalance. As a compromise be- 
tween the superfiuid pairing and an imposed imbalance, 
at intermediate values of the latter, the superconducting 
order parameter condenses at a set of finite center-of- 
mass momenta determined by the details of the Fermi 
surface mismatch and interactions. At sufficiently large 
imbalance, no compromise is possible, and the resonant 
gas transitions to the normal state. 
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FIG. 1: A mean-field zero-temperature phase diagram from 
Refs l79l llOOl of an imbalanced resonant Fermi g a func- 

tion of the inverse scattering length and normalized species 
imbalance P = (A^^ - Ni)/N = AN/N, showing the mag- 
netized (imbalanced) superfluid (SFm), the FFLO state (ap- 
proximated as the simplest FF state) confined to a narrow 
red sliver bounded by Ppflo and Pc2, and the imbalanced 
normal Fermi liquid. 

As illustrated in Figjl] the key observation is that, 
despite strong interactions, within simplest mean-field 
treatments the conventional FFLO state[5ni [ST] remains 
quite fragile, confined to a narrow sliver of polarization 
in the BCS regime[5l |79l \lUO\ HSj. I emphasize that 
in fact above conclusion is only rigorously valid for the 
FF and not other forms (e.g., LO) of the FFLO class of 



states. Although the upper boundary, hc2, just below the 
normal state is trustworthy, as it is shared by all FFLO 
states, the lower one, hd can strongly depend on the form 
of the FFLO state, but was determined by Sheehy and 
Radzihovsky only for the FF stated [1001 EH]- Further- 
more, motivated by earlier studies of the Bogoluibov-de 
Gennes (BdG) equation for the LO state [134Hl36j . com- 
bined with finding of a negative domain- wall energy in an 
otherwise fully-paired singlet BCS superfluid in Zeeman 
field |136| I137j . these studies have quite convincingly ar- 
gued, that a more generic pair-density wave state (that 
includes a larger set of coUinear wave vect or s 1 1 1 8 1 1119p 
may be significantly more stable. The quantitative ex- 




FIG. 2: An illustration of a continuous commensurate- 
incommensurate (CI) transition at hd from a fully- 
gapped (balanced) paired-superfluid to an imbalanced Larkin- 
Ovchinnikov superfluid. The excess majority atoms are local- 
ized on the domain walls in (zeros of) the LO order parameter, 
whose number Udwih) is then proportional to the imbalance 
P{h) and grows continuously with the chemical potential dif- 
ference (Zeeman energy), h — hci- 

tent of the energetic stability of a PDW states in the 
imbalance-detuning phase diagram, in my view remains a 
widely open and urgent question. Consistent with above 
arguments but in absence of controlled quantitative anal- 
ysis, I take the optimistic point of view that the LO states 
can extend over significantly wider region of the phase 
diagram, as schematically illustrated in Fig|3] Assum- 
ing it is indeed energetically stable, its phenomenology 
has been explored beyond its mean- field cartoon [SOI HI] 
(latter only appropriate in the solid state and optical 
lattices, but not in the isotropically-trapped resonant 
atomic gases, where fluctuations are large). 



B. Summary 

Before turning to details, I summarize the salient fea- 
tures of the isotropically-trapped coUinear class of FFLO 
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FIG. 3: A proposed P = AN/N vs. l/(fcFas) phase diagram 
for an imbalanced resonant Fermi gas, showing the more sta- 
ble LO liquid crystal phases (discussed in the text and il- 
lustrated in detail Figjsl replacing a portion of the phase- 
separated regime. 

states, beyond their mean-field approximation. The key 
observation [66l [67] is that such striped states sponta- 
neously break continuous rotational symmetry, and as a 
result exhibit phonon-like Goldstone modes that are qual- 
itatively (energetically) softer than the solid-state analogs 
(where only discrete rotational symmetry can be bro- 
ken) and all other supcrfluids. Namely, striped FF |138| 
and LO classes of states are characterized by highly 
anisotropic (with the modulation wavevector qo = qoz 
and _L transverse to it) collective spectra, with 

uj{k^,k.) - ^jKki + Bkl (1) 

rather than by the usual linear-in-momentum Bogoluibov 
sound mode. The more experimentally relevant striped 
LO state, also exhibits a quantitatively anisotropic Bo- 
goluibov linear-in-fc sound mode, but with the superfluid 
stiffness ratio p']^/ pi that vanishes as the hc2 transition 
to the normal imbalanced Fermi gas is approached from 
below. 

As a result, the fluctuations in such "soft" superfluid 
smectic states are qualitatively stronger. Although the 
states are stable to quantum fluctuations, in 3d the LO 
and FF long-range orders are marginally unstable at any 
nonzero temperature. Consequently, (seemingly para- 
doxically) inside the LO state the average LO order pa- 
rameter 

(Aio(r))fl = (2A,„e'^Wcos(qo-r + 0(r)))«, 

~ ;^cosqo-r — > 0, (2) 

vanishes in the thermodynamic limit (of a large cloud 
with atom number N and cloud size i? — >■ oo), sup- 
pressed to zero by thermal phonon 9/qo fluctuations. 
The LO state is therefore strictly speaking homogeneous 
on long scales, exhibiting "algebraic topological", but no 
long-ranged translational order. Namely, the mean-field 
approximation fails qualitatively and the state instead 
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is characterized by power-law order-parameter correla- 
tions, distinguished from the spatially short-ranged dis- 
ordered phase by confined topological defects (bound dis- 
locations), not by a nonzero LO order parameter. It is 
therefore a 3d analog of the more familiar quasi-long- 
range ordered superfluid film, a 2d easy-plane ferromag- 
net and a 2d crystal [TITifm] . 

As a consequence, a 3d LO state is characterized by a 
static structure function S{q) and momentum distribu- 
tion function n(k) with universal anisotropic gwasi-Bragg 
peaks (around qo), akin to the Landau-Peierls P^Ul 1140] 
behavior of films of a conventional superfluid and 2d 
crystals [l41H144j . Such novel behavior is not, however, 
exhibited by 3d crystalline FFLO states with multiple 
non-coUinear ordering wavevectors jllHl 1119) . that, in 
contrast are characterized by the long-range positional 
order and a nonzero pair-condensate, that is stable to 
thermal fluctuations. 

Another fascinating feature that arises because the LO 
order parameter, A^o is a product (rather than the sum, 
as it is in a conventional supersolid) of the superfluid and 
density- wave component, is the unusual topological exci- 
tation that is a half-vortex bound to a half-dislocation - 
in addition to integer vortices and dislocations. 

In 2d, at nonzero T the LO state is even more 
strongly disordered, characterized by short-range posi- 
tional order with Lorentzian structure function peaks, 
and unstable to proliferation of dislocations 145 . The 
state that results from such dislocated superfluid smec- 
tic is either a "charge" -4 (paired Cooper pairs) nematic 
superfluid 37, 66 or a nematic (possibly "fractionalized" ) 
Fermi liauid |35i 147] . latter qualitatively the same as the 
deformed Fermi surface state [52] . 

Furthermore, a consideration of states that arise due 
to unbinding of various combination of topological de- 
fects (illustrated in the flow-chart in Figj4]) leads to a 
rich array of LO descendent states, that generically must 
intervene between the LO superfluid and a fully-paired 
conventional (isotropic and homogeneous) superfluid and 
a conventional polarized Fermi liquid. If indeed, as ar- 
gued above, the 3d LO state is energetically stable, these 
novel states are expected to appear in the region col- 
lectively denoted "LO liquid crystals" of the detuning- 
polarization phase diagram of FigjSj They include a 
nonsuperfluid smectic {FL?g^_^, driven by an unbinding 
of integer 27r- vortices) , and a superfluid {SFj^, driven by 



a proliferation of integer a-dislocations) and a nonsuper- 
fluid {FLiq, driven by an unbinding of both vortices and 
dislocations) nematics, and the corresponding isotropic 
states, when disclinations also condense. In addition, 
a variety of topologically-ordered isotropic and nematic 
"fractionalized" Fermi-liquid states [FL*j^, FL*^, FL}, 
and others) were predicted [SHI [S3' ^^sA are distinguished 
from their more conventional fully-disordered forms by 
gapped (bound) half-integer defects. These phases are 
summarized by a flowchart Fig(4| and a schematic phase 
diagram illustrated in Fig[5| 

Finally, the ferniionic sector of the LO gapless super- 




FIG. 4: A flowchart of superfluid (SF) and nonsuperfluid 
(FL) phases, exhibiting smectic (Sm) and nematic (TV) con- 
ventional orders as well as topological orders (indicated by * 
and **), induced by a proliferation of various combination of 
topological defects, (0,a), (27r,0), and (n,±a/2). 



conductor is also quite unique, exhibiting a Fermi sur- 
face of Bogoliubov quasiparticles associated with the An- 
dreev band of states, localized on the array of the LO 
domain walls. Consequences of the interplay between 
these fermionic and Goldstone mode degrees of freedom 
remains an open problem. 



C. Microscopic model of imbalanced Fermi gas 



grand-canonical Hamiltonian 



The physics of imbalanced atomic Fermi gases interact- 
ing through a broad Feshbach resonance is well captured 
by the one-channel model [SHSl I146| , characterized by a 



k,(T 



A<rr)cLc, 



ko-^ko- 



kqp 



with the single-particle energy — h^k^ /2m. The sepa- 
rately conserved number N^r — {N^, N^) of atomic species 
(hyperfine states) a ='f , i) are imposed by two chemical 
potentials, fj-^r = {p.-f,fJ.i). 
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FIG. 5: A schematic imbalance-chemical potential (Zeeman energy), h — jj.-^ ~ /i^ ys detuning (interaction strength), —l/kpa 
phase diagram, illustrating the 3d LO smectic phase (SFsm) and its descendant (described in the text), driven by a proliferation 
of various combinations of topological defects. The inset shows the global imbalance-interaction BCS-BEC phase diagram, 
illustrating the location of these putative phases. 



The key feature that distinguishes this Fermi system 
from those familiar from solid state contexts is the attrac- 
tive resonant interaction parameterized by a short-range 
s-wave pseudopotential, g < 0. Through an exact T- 
matrix scattering calculation g controls the magnetic- 
field tunable [148) scattering length 
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that diverges above a critical attraction strength, g^ — 
(A ^ 1/d is the short-scale pseudo-potential cutoff 
set by the extent of the molecular bound state), corre- 
sponding to a formation of a two-atom bound state. 

The many-body thermodynamics of the resonant 
Fermi g function of N,AN,T,as, (i.e., the ex- 

tension of the BEC-BCS crossover to a finite imbalance, 
AN) at large /cFlOisI presents a formidable challenge. 
However, much progress has been made in mapping out 
its qualitative (and in some regimes quantitative) phe- 
nomenology through a variety of approximate theoretical 
techniques, including quantum Monte Carlo "TW , mean- 



field theory [Ti HOI HOOl HOD [M] , the large- A^/ (fermion 
fiavor) [inniinn] and e-expansions [W] . 

The simplest of these is the standard mean-field 
analysis [7^ llOOj that gives a satisfactory qualitative de- 
scription (quantitatively valid deep in the weakly-coupled 
BCS regime, fcirjasl <C 1), as a starting point of more 
sophisticated treatments. To this end we assume the ex- 
istence of a condensate 



A(r) = En 



3«q r 



g(c;(r)c^(r)). 



(5) 



corresponding to pair-condensation at momenta q, with 
the set of amplitudes Aq and q to be self-consistently 
determined by the minimizing the ground state energy 
subject to the constraints of fixed total atom number 
N = + and the atom species number imbalance 
("polarization") AA^ — — N^, imposed by the aver- 
age and difference chemical potentials n,h — |(/i^ ± ni), 
latter corresponding to the pseudo Zeeman energy. 

Specializing to the simplest case of a single q of the 
Fulde-Ferrell state [SUI this reduces the Hamiltonian to a 



6 



quadratic Bogoluibov form, that can be easily diagonal- 
ized. This gives [79 t llOOj the ground state energy 



takes a famihar form 



and the excitation spectrum i^kcr 
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and El- 



(7) 



Eqq (and its generahzations to finite-temperature free 
energy [TD3]), then gives all the thermodynamics, includ- 
ing the phase behavior summarized by the phase diagram 
in Figjl] 

In particular, this analysis predicts the existence of the 
FF state, stable only over a narrow sliver of imbalance, 
closing down for ~l/{kFa) > 0.5 [79l llOOj . As men- 
tioned in the Introduction, there are compelling argu- 
ments suggesting that this fragility is specific to the single 
q planewave FF condensate, and that the more generic 
PDW states are far more stable because they allow ener- 
getically important amplitude modulation [HZl ll34H137j . 
However, general PDW states are difficult to analyze at 
the transition from the fully gapped, balanced paired 
superfluid, near the lower critical field /id (at vanish- 
ing imbalance). In contrast, an analytic analysis near 
the upper-critical chemical potential difference hc2 at the 
transition from the normal state is indeed possible as a 
controlled Ginzburg-Landau expansion in the small pair- 
density wave amplitude, Aq [SH [Ml ISZl [HZ] • While not 
quantitatively accurate away far below the hc2 transition 
(where PDW order parameter is large and is not limited 
to a single Fourier component A^) such Landau expan- 
sion is expected to be qualitatively correct and is a good 
starting point for a more complete analysis of fiuctuations 
and phase transitions into the PDW LO state. 



D. Order- parameter theory of FFLO states 

1. Ginzburg-Landau expansion near hc2 

The analytical treatment of the FFLO states near hc2 
relies on the Ginzburg-Landau expansion in Aq, that is 
small near the (in mean-field) continuous hc2 normal-to- 
FFLO transition [5T1 1149] . This expectation is supported 
by the exact Id BdG solution[134] at high fields, where 
A(a;) is indeed well-approximated by a single sinusoid, 
with an amplitude A^ that vanishes continuously near 

hc2- 

Consistent with these general arguments, by integrat- 
ing out the atomic degrees of freedom, near hc2 the 
Ginzburg-Landau expansion for the ground-state energy 



q 



:,|Aq| 



{q.} 



A* A A* A 
q2 ,q3 ,q4 ^qi "q2 ^qa ^q4 ; 



(8) 



(6) where the dispersion is given by[5Tl [67l [79l llOOj 



3n 



-1 



^BCS 



4h^ 



In 



2h 



(9) 



whose minimum at a finite qoih) w 1.2h/vp (near hc2) 
captures the imbalanced atomic Fermi system's energetic 
tendency to pair at a finite momentum, and thereby to 
form a pair-density wave characterized by a reciprocal 
lattice vector with magnitude qg and a spontaneously 
chosen orientation. The value of h at which Sq^ van- 
ishes determines the corresponding mean-field N-FFLO 
transition point. While at quadratic order, all Fourier 
modes with |q| = qo are degenerate, becoming unsta- 
ble simultaneously, the form of the FFLO state is dic- 
tated by the interaction vertex function, Vq^.qj^q^^q^ that 
has been explicitly computed. Near the transition the 
physics of a unidirectional pair-density wave (Cooper- 
pair stripe) order, characterized by a collinear set of q„'s 
is well captured by focusing on long-wavelength fiuctua- 
tions of these most unstable modes, well described by a 
Ginzburg-Landau Hamiltonian density 

n = J [iV^Ap - 2q2|VAp] +r|A|2 + h,\A\^ + h^i^ 

(10) 

where deep in the BCS limit, near the hc2 the model 
parameters are given by 



J 



OMn 

— ^1 qo 

epqo 



-A 



BCS, 



Ai 



l.SlAacs 
3n 



3n 
4ep 



In 



, (11a) 



BCS 



9h 

4:hc2 

1-83W ^^^^^ 
^F^licslo ' 



and the inclusion of the current-current interaction, 
j = ;^Re [— A*(r)zVA(r)] is necessary for a complete de- 
scription. More generally, away from the weak-coupling 
BCS limit these couplings can be taken as phenomenolog- 
ical parameters to be determined experimentally, but the 
general form of the Ginzburg-Landau model has broader 
range of applicability, capturing all the qualitative fea- 
tures of the transition and the PDW state. 



2. Larkin-Ovchinnikov state near hd 

However, the derivation and expressions for the associ- 
ated couplings, Eqs.( 11 1 are limited to the weak coupling 
BCS regime and near the high chemical potential imbal- 
ance (Zeeman field) normal-to-FFLO transition at hc2- 

In a complementary, low chemical potential imbalance 
regime, just above the transition from the fully-paired 
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(BCS-BEC) superfluid to the LO state at hci, a phe- 
nomenological analysis is possible[67]. It treats the LO 
state as a periodic array of fluctuating ±A domain-wafls 
(stripes) in A(r), akin to the lyotropic phases in soft con- 
densed matter [48, i49j . 

However, such approach imphcitly assumes that as the 
domain- wall surface energy becomes negative [T31HI3Z] 
for h > hci, their interaction remains repulsive, and so 
the domain-walls proliferate continuously as a periodic 
array inside the LO state. Under this assumption (that 
warrants further study) the domain- wall density n^to and 
the associated species imbalance P oc Udw (~ <lo{h)) 
is then set by a balance between the negative surface 
energy and the domain-wall repulsion, growing continu- 
ously as a function of h — hd according to the Pokrovsky- 
Talapov's commensurate-incommensurate (CI) transi- 

This behavior is clearly ex- 



E„(A, h) 



tion phenomenology [150] 
hibited in ld jl341 I152[ I153 j through an exact solution 
and bosonization methods, and has been argued to per- 
sist in higher dimensions |134H137] . The CI route for 
a transition to the LO state contrasts sharply with the 
Landau theory [STl 17^ 1100] of two independent order pa- 
rameters Aq, Ag, that always predicts a first-order BCS- 
LO transition. The latter corresponds to the case of an 
attractive domain-wall interaction, that therefore prolif- 
erate discontinuously above hd , leading to the ubiquitous 
phase separation found in mean-field theorv [7^ 1100] . It 
is currently unclear what dimensionless microscopic pa- 
rameter, analogous to Abrikosov's k (distinguishing be- 
tween type I and type II superconductors) 154, 155,,, con- 
trols these two alternatives of the macro-phase separation 
(a first-order transition) and the micro-phase separated 
LO state (a continuous transition out of the gapped SF 
state) [133] ■ A detailed analysis of such low-imbalance 
approach to the SF-LO transition and the LO state is 
sorely missing and is a subject of current research. 

A semi-phenomenological local density approximation 
(LDA) model that assembles all known ingredients is 
given by 



H[Aiv)] 



IV^AP 



2g2|VA|2) + y(A(r)) 



(12) 



where Aq = c^/'^Abcs and V{A) = -2v{^i)A'^ In (^) - 

v{^jl) [hy^h^ - A2-A2 cosh"i(VA)] Q{h-A) is the effec- 
tive potential derived within a BdG analysis for a uniform 
A. It fully captures the double-minimum structure and 
the associated Ist-order normal to (fully-gapped BCS) 
superfluid transition that skips the interesting interme- 
diate e.g., the FFLO states. This LDA potential is sup- 
plemented by the gradient energy functional, inherited 
from a microscopic GL analysis (valid only near /ic2)- Its 
use near hci is supported by the fact that it is the sim- 
plest form that incorporates the underlying symmetries 
and encodes the expected energetics for the system to 
order at a finite /i-dependent momentum even near hd- 
Work is under way to derive this functional through a 




A 
Abcs 



FIG. 6; Ground state energy Egs[A,h] as a function of the 
order parameter A and Zeeman energy h, indicating a first- 
order transition at he (thick red curve) between the paired su- 
perfluid state (Abcs) and the normal state (A — 0), at fixed 
density and/or imbalance exhibiting phase-separated coexis- 
tence. As argued in the text, analogously to the Pokrovsky- 
Talapov systems and type-II superconductors, this first-order 
transition can be preempted by a continuous Cl-like transi- 
tion to a striped LO superfluid and its descendent quantum 
liquid-crystal states, illustrated in Figj2] 



controlled Moyal (semi-classical) gradient expansion on 
the BdG Hamiltonian. 

The functional EQ[A{r)] has a double- well structure, 
with an additional normal state minimum (at A = 0) de- 
veloping for h > ABC's/2. It thus allows periodic soliton 
structure in A(r), corresponding to oscillations between 
the minima at ±Abcs- 



E. Goldstone modes in striped FFLO states 

Using the Ginzburg-Landau model one can develop a 
low-energy Goldstone modes description of the striped 
paired states and use it to analyze their stability to fluc- 
tuations. We focus on the unidirectional (striped) pair- 
density wave states, with FF and LO states as simplest 
representatives of two qualitatively distinct universality 
classes. The corresponding order parameter is given by 



AfFLo(r) = A+(r)e^i- + A_ 



.(r)e' 



zq-r 



(13) 



where A± (r) are two complex scalar order parameters, 
the dominant Fourier coefficients A-|-(r) = A^(r)e*'^='=('') 
and amplitudes Aj. distinguishing between the FF and 
LO states. Using this representation inside H and min- 
imizing the ground-state energy for h < hc2 a sim- 
ple analysis shows that indeed it is the LO state with 
A+ = A_ 7^ that is most stable inside the BCS regime; 
the FF state is characterized by only one nonvanishing 
order parameter, A+ ^ 0, A_ = 0(51]. More generally, 
either state can be stabilized depending on the relative 
magnitudes of Ai and A2. 
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1. Fulde-Ferrell state 

The FF state is characterized by a single (independent) 
nonzero complex order parameter, 

AFp{r) = A,,e'^o-r+i^, (14) 

that is a plane-wave with the momentum qo and 
a single Goldstone mode = 0+ corresponding to 
the local superconducting phase. The state carries a 
nonzero, uniform spontaneously-directed supercurrent 
jpp = i|Agop(qo -I- V0) and thereby breaks the time- 
reversal and rotational symmetry, chosen spontaneously 
along qo, as well as the global gauge symmetry, corre- 
sponding to the total atom conservation. Although the 
FF order parameter itself is not translationally invari- 
ant, it is invariant under a modified transformation of 
an arbitrary translation followed by a gauge transfor- 
mation, with all gauge-invariant observables thus trans- 
lationally invariant. Thus, the FF state is a uniform 
oricntationally-ordcrcd polar superfluid. The underly- 
ing rotational invariance also demands that it is in- 
variant under a rotation of qo = qoz by an angle a 
that generates a nontrivial, spatially-dependent phase 
(p'^ir) = z(cosa — 1) -I- xsina. Simple algebra demon- 
strates that the fully nonlinear form of the longitudinal 
current = d\\4' + ^Qo^i'^'t')'^ ensures that it and the 
corresponding energy HpF vanish for 0'^(r), as required 
by the underlying rotational invariance. 

The analysis of the GL functional lead to a Goldstone 
mode (j) Hamiltonian|66 1 \67 \ I138j 

-Hff = lKiV^cbr + lpl{d\\c^+lqo\y<Pr)\ (15) 

where 9j| = qo • V, pL' ~ 8Jq^|Aq„p is the superfluid 
stiffness along qq and K ~ 2J\Aq^\. The Hamiltonian 
form, Hff is valid beyond its weak-coupling microscopic 
derivation and is familiar from studies of conventional 
smectic hquid crystals [48 l l49 l [T56] . despite the fact that 
FF state is a translationally-invariant polar superfluid 
not a smectic. The necessity of keeping the higher or- 
der gradients and nonlinearities inUpp is due to the 
identical vanishing of the transverse superfluid stiffness, 
pj- = (guaranteed by the underlying rotational invari- 
ance unique to the FFLO striped states, absent in solid 
state contexts) that leads to fluctuations that are other- 
wise infrared-divergent in a purely harmonic model. 

8. Larkin-Ovchinnikov state 

The LO state is instead described by two independent 
nonzero PDW amplitudes A-|_ = A_ = A^j, (growing 
below hc2), that lead to a standing wave pair-density 
wave order parameter, 

Aio(r) - 2A,„e*5W++*-)cos[qo-r-Fi(0+-0_)],(16a) 
= 2A,„e*^cos[qo-r + 0], (16b) 



that is a product of a superfluid and a unidirectional 
density wave striped order parameters. These are re- 
spectively characterized by two Goldstone modes cf),6, 
corresponding to the superfluid phase and the smectic 
phonon u = —O/qq of the striped state. This also con- 
trasts with the conventional smectic[48, (e.g., in liquid 
crystal materials, where one instead is dealing with a real 
mass density p(r) not a pair condensate wavefunction) , 
characterized by a single phonon Goldstone mode, u. 

The mean-field LO order parameter, A^o simultane- 
ously exhibits the ODLRO (superfluid) and the smec- 
tic (unidirectional density wave) orders. It thus sponta- 
neously breaks the rotational, translational, and global 
gauge symmetries, and therefore realizes a form of a 
paired supersolid. However, it is distinguished from a 
conventional purely bosonic supersolid [125H128"] . where 
homogeneous superfluid order and periodic density wave 
coexist, by the vanishing of the ("charge"-2 two-atom) 
zero momentum (q = 0) superfluid component in the LO 
condensate jl29j . 

Similarly to the FF state, the underlying rotational 
symmetry of the LO state strongly restricts the form of 
the Goldstone-mode Hamiltonian. Namely, its 9 = —qqu 
sector must be invariant under a rotation of qo, that 
defines the spontaneously-chosen orientation of the pair- 
density wave, and therefore must be described by a smec- 
tic form [151 I156j . On the other hand because a rota- 
tion of the LO state leaves the superconducting phase, 4> 
unchanged, the superfluid phase sector of the Hamil- 
tonian is therefore expected to be of a conventional xy- 
model type. Consistent with these symmetry-based ex- 
pectations the LO Goldstone-mode Hamiltonian was in- 
deed found [SSJ in?] to be given by 

-Hlo = iif(VV + ^B(9p-i(Vt.)2)' 

+ ^P^(9||0)' + ^P^(V^0)2, (17) 

with the nonlinear strain tensor Uqq — q • Vu — ^(Vu)'^ 
ensuring the full rotational invariance. In the weakly- 
coupled BCS limit the smectic elastic moduli and the 
superfluid stiffnesses are given by 

K = 4Jql\Aqf ^ ^■^''^Bcs i^^h/h,,), (18a) 

fF^o 

B ^ 16Jgo"|A,„p^ ^-^"^'^^^ ln(fe/fec2), (18b) 
pI = B/ql (18c) 

Pj = ^-^lAqX-^'^'^^^Ah/K,). (I8d) 

Thus, the LO state is a highly anisotropic superfluid 
(though less so than the FF state, where =0), with 
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a ratio that vanishes for h — >■ h~2 ■ 

We stress that while the detailed expressions for the 
moduli above are specific to the weak-coupling BCS limit 
near hc2 the general form of Hlo, ( |17| ), including the 
structure of the symmetry-enforced nonlinearities in the 
u (smectic) sector is valid beyond our microscopic deriva- 
tion, and holds throughout the LO phase. 

By extending the Hamiltonian to include density fluc- 
tuations, Sn±, canonically conjugate to 4>± and integrat- 
ing them out in an imaginary time (r) coherent-state 
path integral, leads to a Lagrangian density 



C 



Xo 



(20) 



with H given by the Hlo in the LO ground state (FF 
state is treated similarly using HpF and a single Gold- 
stone mode). For the LO state, this analysis then 
predicts the existence of two anisotropic low-frequency 
modes with dispersions 



W0(k) = \/ {p^kl+ plk^J/xo, 



(21a) 
(21b) 



where xo is the compressibility of the Fermi gas. These 
modes respectively correspond to the zeroth sound (the 
Bogoliubov mode as in a conventional superfluid) and 
smectic phonon, unique to the LO state. In cold atomic 
gases, these should in principle be measurable via the 
Bragg spectroscopy techmaue |157H15iJ| 

With the Goldstone-mode Lagrangian in hand, the 
effects of quantum and thermal fluctuations as well as 
equilibrium correlation and response functions can be 
calculated [TBD]. 



3. Goldstone modes fluctuations 

Armed with the action for the Goldstone modes, the 
low-energy quantum and thermal fluctuations in the 
FFLO striped states are straightforwardly computed. 
Despite smectic like softness of these modes, it is easy to 
show that quantum fluctuations remain finite for d > I 
and therefore the FF and LO states remain stable at zero 
temperature. 

Harmonic approximation: 

In contrast, at finite temperature the root-mean- 
squared fluctuations of the smectic phonon modes in the 
LO state (and phase fluctuations in the FF state) diverge 
with trap size (in a purely harmonic Goldstone-modes 
theory) according to 



T r3-rf 



4ttVBK 



d<3, 
d = 3, 



(22) 



where the fluctuations are evaluated in a trap with an 
aspect ratio L± x L^, with its extent along the or- 
dering qo axis [IHl EHl I162j . These fluctuations lead 



to an emerge nce of important crossover length scales 
^z,C± ^ i^z\/K/B)^^'^ = "v^^^A that characterize the 
finite-temperature LO state. 



/BK 



T 



K 

^ Tqo' 
BK/T , 



d 

cK 



2, 
3, 



(23) 



defined as scales at which LO phonon fluctuations are 
comparable to LO stripe period a — 2Tr/qQ. 

The thermal connected correlation function of LO 
phonons 



(24) 



is also straightforwardly worked out, in 3d giving the 
logarithmic Caille form [162] 



Cf(r^,z) 



T 



27rVKB \ In(^), 




(25) 



In 2d it is instead given by 



Cf(x,z) 




X < 
X > 




(26) 



Above finding of the divergence of smectic phonon fiuc- 
tuations at nonzero temperature have immediate drastic 
implications for the properties of the LO (and FF) states. 
The most important of these is that the thermal average 
of the Landau's LO order parameters (16) vanishes in 
thermodynamic limit 



(Aio(r))o « 2Ag„(ii)cos(qo-r), 



(27) 



with the thermally suppressed order parameter ampli- 
tude given by 



n/2 



d = 2, 
d = 3, 



-> 0, for Lj_ 



(28) 



where 



accounts for the finite quantum 



and thermal superconducting phase fluctuations, and 



V 



is the Caille exponent |162j . 



SwVBK 

Akin to 2d xy-model systems, the vanishing of the LO 
order parameter does not imply the instability of the 
phase, but that the true fluctuating state contrasts qual- 
itatively with perfectly periodically-ordered mean-field 
description. 

These divergent LO phonon fluctuations also qual- 
itatively modify the Cooper-pair momentum distribu- 
tion function rik = (Aj^Ak) and structure function S'(q) 
(respectively measurable via time-of-flight and Bragg 
spectroscopv |157Ul59j imaging of pair-condensate) from 
their true Bragg (5-function peaks) form characteristic of 
the mean-field long-range periodic order. They are highly 
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LO order is even more strongly suppressed by thermal 
fluctuations down to short-ranged correlations, strongly 
suggesting true instability of the LO order. Because it 
is the "soft" smectic Goldstone mode that is responsible 
for these interesting properties they are necessarily also 
shared by the FF state[n3- 

I emphasize that analogous Kosterlitz-Thouless phase 
fluctuation physics has been observed in conventional 2d 
trapped superfluids, despite the finite trap size |164lll65j . 
Thus, I am hopeful that they can be similarly seen in the 
3d LO state. A more detailed analysis of the trap effects 
is necessary for direct comparison with experiments. 

Another fascinating consequence of the vanishing 
of (Alo)i (28) is that the leading nonzero Landau 
order parameter characterizing the LO state is the 
translationally-invariant "charge" -4 (4-atom pairing) su- 
perconducting order parameter. 



FIG. 7: The finite momentum pairing at qo and divergent 3d 
smectic phonon fluctuations in the LO state are reflected in 
the Cooper-pair center-of-mass momentum distribution func- 
tion, rik (observable via the time-of-flight measurements) , dis- 
playing power-law Bragg peaks, characteristic of spatial quasi- 
long-range order. 




2qo 4qo 6q, 







FIG. 8: The structure function, S'(q) for the 3d LO state, 
displaying power-law Bragg peaks, characteristic of the LO 
superfluid's spatial quasi-long-range order. 



anisotropic {q^ ^ ^i/A) and exhibit quasi-Bragg peaks 
(see Figjs]) around the ordering wavevector qo (and its 
harmonics, q„), reminiscent of (l-l-l)d Luttinger liquids 
and two-dimensional crvstals [139H144j 



LO 



(29a) 
(29b) 



These predictions are a reflection of the well-known [48l 
and experimentally tested [163| behavior of conven- 
tional smectic liquid crystals. In two dimensions, the 



A(4) 



(c,(r)ct(r)c,(r')ct(r')) 



/\\ 'iq-(r-r') 



r,r' 



-kCt,k+qCi -k'Cf.k'-q) 

k,k' 



(30) 



Thus in the presence of thermal fluctuations the LO 
phase corresponds to an exotic state in which the off- 
diagonal order is exhibited by pairs of Cooper pairs, i.e., 
a bound quartet of atoms, rather than by the conven- 
tional 2-atom Cooper pairs [371 ES]- In 2d and 3d this 
higher order pairing is driven by arbitrarily low-T fluctu- 
ation, rather than by a fine-tuned attractive interaction 
between Cooper pairs, and therefore has no simple mean- 
field description. A microscopic formulation of such a 
state and its detailed properties remain an open prob- 
lem. 

Nonlinear elasticity beyond Gaussian fluctuations: 
As discussed in the context of conventional 
smectics |156j and more recently for the FFLO 
state[Bni EZ], above predictions neglect the ef- 
fects of Goldstone mode nonlinearities 'Hnoniinear = 
-iB(a^u)(VM)2 + |B(Vu)^ in Eq.((l5|,([T7|, that mod- 



ify the asymptotics on scales longer than the crossover 
scales = a/CF^, 



NL 



1 / Kl 

T \ B 



1/2 



, d = 2, 
, d = 3, 



(31) 



The behavior on scales beyond ^^^^ can be obtained us- 
ing renormalization-group analysis for d < 3 |671 I156j . 
with an exact solution in 2df l66j . The finite-temperature 
asymptotics is well-approximated by a correlation func- 
tion 



G„(k) 



T 



B(k)fc2+X(k)fc4 ^ 



(32) 



with moduli B(k) and if (k) that display a universal sin- 
gular wavevector-dependence, that is asymptotically ex- 
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act logarithmic [TS3] in 3d 



55 



K{k^,k,^0) ^ i^|l + -^ln(l/fcia)p/5 , (33a) 

D47r 



B(kj_ = 0,k,) - B\l 



and power-law in 2d 



5g 

1287r 



ln(A/fc,a^)| 



2n|-4/5 



(33b) 



K{k) = if(fc^eD , (34a) 



B(k) ^ B [k^e^^y^ fB{k.s.ri{k^irf) , (34b) 

with fsix), fxix) universal scaling functions and rj^g = 
1/2, 77!^ = 1/2, C^'^ = 3/2 exact [US]. In 3d this translates 
into an equal-time LO order parameter correlations given 
by[ISS] 

n(z,ri = 0) = (A^o(z)Aio(O)), (35a) 
- e-'=i('"^)°^'cos(go2). (35b) 

Although these 3d anomalous effects are less dramatic 
and likely to be difhcult to observe in practice, theo- 
retically they are quite significant as they represent a 
qualitative breakdown of the mean-field and harmonic 
descriptions of the FFLO striped states. 



the LO state admits two types of topological defects, 
characterized by integers 7V„ = (?T._|_,n_) defined by 
§ d£ ■ V(/)± = 27rn±. These equivalently correspond to 
superfluid vortices and edge dislocations in the striped 
PDWfTBT]. These are characterized by multiples of 
/ia//-integers ny,nd = (n_|_ ± ri_ ) /2 and therefore allow 
four types of elementary defects: integer (27r, 0) vor- 
tex, integer dislocation (0, a) and two half-integer vortex- 
dislocation composites (tt, ±a/2). The latter composite 
fractional defects are allowed because a sign change in 
A^o due to a a/2-dislocation in u is compensated by a 
TT-vortex in (j) (thereby preserving a single-valuedness of 
ALo)E3liniEIlEZl- In terms of the two coupled , 0- 
Goldstone modes, these correspond to an integer vortex 
in one and no vortex in the other superfluid phase. 

Their thermodynamics and correlations can be treated 
via a mapping on a multi-component Coulomb gas 



jjxy—sm 
^CG 



piPs 

^ q 



''q. t; I 



"T-o d\ 

q 



where 



q 

-psm 

q 



\J pi I pi + ql\J p\/p. 



q\ 



_L 

s 1 



(36) 



(37a) 
(37b) 



Phases and transitions 



Topological defects 



Associated with its two compact Goldstone modes, 0± 
(equivalently 0(r) — |(0+ +</'-), u{v) — \a{4>+ — 4>-)) 



with TOq,^, m^^d the Fourier transforms of the vortex and 
dislocation densities. Equivalently, defects thermody- 
namics can be analyzed via duality transformation. In 
two dimensions it leads to a sine-Gordon-like model for 
dual 0, fields characterizing fractional defects, with the 
dual IIamiltonian,37. ,67. ,171j 



H: 



SG 



i=l0qP 

pi pi 



Kql 



^ql 



ff7r,a/2 C0S( 



COS 7r0 coslTT' 



0) + 5277,0 cos(27r(^) -I- 5o,a cos(27r6i) , (38) 



r 



It is convenient for analyzing the effects of defects on the 
LO state, particularly for a computation of their screen- 
ing on long scales, unbinding, and for the analysis of the 
resulting disordered state. From the form (38) it is clear 



that (aside from an inconsequential anisotropy) the dual 
vortex sector described by (p has a standard sine-Gordon 



form. In contrast, the dual dislocation sector, described 
by 9 is qualitatively modified by the highly nonlocal and 
qualitatively anisotropic smectic kernel, T"^ . 

A standard analysis gives the relative energetics of 
these defects, in the thermodynamic limit (-/j_l,2 ~^ oo) 
given by 



KL<^E 



'(7r,a/2) 



-L\nL 



K 



L < Ey. 



'(27r,0) 



' psLhiL, for L_L ^ ~ L 



oo. 



(39) 
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FLsn. FLn FLi 

Tu(i) Tu(i) Tu(i) 

SFsm ^ SFn ^ SFi 



TABLE I: Five phases that naturally emerge as disordered 
descendants of the LO (superfluid smectic, SFsm) state. 



Based on this energetics one may be tempted to conclude 
that in this limit (unless preempted by a first-order tran- 
sition) it is the integer dislocation loop defects that pro- 
liferate first and the LO smectic preferentially disorders 
into a nematic superfluid, SFj\j. However, in contrast to 
the 2d KT mechanism |144j . the 3d disordering transitions 
take place when the relevant stiffness, renormalized by 
quantum and thermal fluctuations is continuously driven 
to zero at the transition, or takes place at a finite (rather 
than a vanishing) defects fugacity. For a thermal transi- 
tion this roughly corresponds to a transition temperature 

set by the corresponding stiffnesses, Ps = \ff^^ ^^^d 
K, B. Thus, in principle by tuning these stiffnesses via 
imbalance and resonant interaction, a variety of phases 
can be accessed. 



2. Conventional phases and transitions 

By considering all possible basic combinations of spon- 
taneously "broken" subset of spatial and gauge sym- 
metries leads to an array of partially spatially-ordered 
paired superfluids and Fermi-liquid states, that are de- 
scendants of the smectic LO (SFgm) state. These 
isotropic (I), nematic (N) and smectic (Sm) SF and FL 
states are summarized in Table I. ft is notable that the 
isotropic superfluid, SFj exhibits a flnite species imbal- 
ance and off-diagonal long-range order, symmetry-wise 
isomorphic to the polarized superfluid, 5'f at f79 [ll00j . lat- 
ter confined to the BEC side of the BCS-BEC crossover. 
In contrast, (as a descendant of the LO state expected 
to be stabilized by Fermi surfaces imbalance) the SFi 
state is realized in the BCS regime, something that has 
been searched for dating back to Sarma [TT7] . but has not 
been possible within mean-field treatments, that instead 
predict an instability to phase separation|711 [TSl llOOj . 
The isotropic Fermi liquid, FLj is isomorphic to the 
conventional normal state. Together these intermediate 
fluctuation-induced phases naturally interpolate between 
the fully gapped singlet (homogeneously and isotropic) 
BCS superconductor at zero imbalance and low temper- 
ature, and a polarized Fermi liquid at large imbalance 
and/or high temperature. 



phase /symmetry 


U{1) 




R 


FLi 


V 


V 


V 


FLn 


V 


V 


X 


FLsm 


V 


X 


X 


SFi 


X 


V 


V 


SFn 


X 


V 


X 


SFsm 


X 


X 


X 



TABLE II: A summary of LO liquid crystal Fermi-liquid (FL) 
and superfluid (SF) phases, and corresponding order param- 
eters and broken symmetries, indicated by X's. Unbroken 
symmetries (gauge (7(1), translational Tqg, rotational R) are 
marked by check marks. The subscripts I, N, Sm respectively 
indicate the Isotropic, Nematic and Smectic orders. 



3. Topological phases via defects unbinding 

The phases discussed above can be complementarily 
characterized through unbinding of different combina- 
tions of topological defects. The smectic (whether SF LO 
state or FL smectic) to nematic transition is driven by un- 
binding of defects with edge dislocation charge, followed 
by transition into the isotropic state driven by prolifer- 
ation of disclinations. The superfluid and Fermi-liquid 
version of these liquid crystal states are distinguished by 
unbinding defects with superfluid vortex charge. 

However, a characterization in terms of topological de- 
fects also allows a distinction between topologically dis- 
tinct phases with the same conventional order, where 
a Landau order parameter is insufficient to distinguish 
them. In fact because of the vanishing LO order param- 
eter, description in terms of topological order is neces- 
sary even for the smectic LO state, distinguished from its 
more disordered descendants by the absence of unbound 
topological defects, in direct analogy with the quasi-long- 
range ordered state of the 2d xy model. 

A rich variety of possible phases and transitions is dis- 
played in a schematic imbalance-detuning P — 1/kpa 
phase diagram, Fig[5] Increasing the imbalance sup- 
presses the superfluid stiffness and drives the system to- 
ward a conventional Fermi liquid state, FLj at hc2- Con- 
versely, a reduction in species imbalance primarily re- 
duces the elastic moduli of the smectic pair-density wave 
by increasing its period l/go and thereby weakening the 
interaction between the LO domain-walls, driving the 
system toward a conventional isotropic and homogeneous 
superfluid SFj at hd- 

Thus, starting with the LO SFsm state and decreasing 
h leads to the unbinding of the integer dislocations (0, a), 
and a transition to an orientationally-ordered, i.e., a ne- 
matic "charge" -4 superfluid, SF^. The later "charge" -4 
feature of SFj^ naturally appears as the remaining sec- 
ondary order parameter Ai^J = A|^q once the LO po- 
sitional order A^o is destroyed by unbinding of integer 
dislocation loops. Since in contrast, the "charge"-2 SF 
nematic order vanishes in the LO state, a direct transi- 
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tion to it from the LO state can generically only proceed 
through a first-order transition. 

Conversely, increasing h starting with the LO SFsm 
is expected to lead to a suppression of ps, a prolifer- 
ation (27r,0) vortices, and a transition to a 2g-smectic 
Fermi liquid, FL^g^, a non-superfluid periodic state with 
a wavevector that is twice the LO state. Alternatively, 



as suggested by the energetics in Eq. (39), if instead the 



lower-energy half-vortex dislocation defects (tt, a/2) (or 
the (tt, —a/2), but not both) unbinds first, a transition to 
a nematic Fermi liquid, FL*j^ (with the restored transla- 
tional and U(l) charge symmetries) will take place. The 
resulting state is qualitatively distinct from the more con- 
ventional nematic (orientationally ordered) FLj^ phase 
in which both (7r,a/2) and (tt, — a/2) are proliferated. 
Both are also distinct from the nematically ordered FL*j^ 
state, in which only integer dislocations, (0, a) and inte- 
ger vortices, (27r, 0) are unbound. One can envision a 
number of other states and phase transitions at low h 
by further considering the disordering of the nematic su- 
perfluid, SF^^ by unbinding various patterns of disclina- 
tions and 7r-vortices. Many open questions remain about 
the relative energetics and detailed properties of these 
phases. 

Above considerations lead to at least three topologi- 
cally distinct Fermi liquid phases that naturally emerge 
from disordering of the LO (S'-Fsm) phase by unbind- 
ing different combinations of allowed defects. Because 
the conventional vortex (27r, 0) and the conventional dis- 
location (0, a) are composites of the fundamental defects 
(tt, ±a/2), the nonsuperfiuid states FL*j^, FL*^ and their 
isotropic cousins FL*j, FL" (in which disclinations are 
also unbound) are expected to be "fractionaUzed" [T75] . 
topologically distinct from their conventional Fermi liq- 
uid analogs, where (7r,±a/2) are also unbound. 

These novel phases are analogous to the putative 
phase-disordered fractionalized states obtained by un- 
binding double (hc/e) vortices, studied extensively 
by Sachdev, and by Balents, Senthil, Fisher, and 
collaborators 1 1 72T11 74j in the context of high tempera- 
ture superconductors. The resulting nonsuperfiuid phase 
is distinguished from a conventional Fermi liquid by a 
gapped "vison" , a Z2 defect that is a remnant of the fun- 
damental hc/2e vortex after the composite hc/e (double) 
vortices proliferate. 

The states FL*j^, FL*^ also bare a close relation to the 
collective mode fractionalization discussed by Sachdev, 
et al. |175fll77] in the context of quantum paramag- 
netic phases, emerging from disordering a collinear spin- 
density wave. As with the (C/(l) ® U{1)) /Z2 LO state, 
where the order parameter is a product of the superfluid 



and smectic order parameters, Eq.(16l, there too the or- 



der parameter is of (6*2 (81 C^(l)) /Z2 product form, en- 
coding spatial modulation of the spin density, and there- 
fore admits half-integer (vison-like) defects. Correspond- 
ingly, the phases and transitions can equivalently be cap- 
tured via an effective Ising gauge theory. Its nonzero Z2 
flux through a plaquette encodes the presence of a half- 



integer (tt, a/2) defect. The Ising gauge field encodes the 
loca l Z2 redundancy of splitting the LO order parame- 
ter (16) into a "charge"-2 boson, 6j — e^"^', that cre- 



ates a zero-momentum Cooper-pair (diatomic molecule) 



and a neutral boson, p| 



t 



that creates a den- 



sity wave at the LO wavevector q. In the above notation 
this is the nonsuperfiuid periodic state dubbed FL^g^, 
in which (27r, 0) vortices have proliferated, but disloca- 
tions remain bound. For a vison remaining gapped, the 
resulting nonsuperfiuid nematic state is the topologically 
ordered FL*j^, qualitatively distinct from a conventional 
FLjv in which the vison is gapless. The two FL phases 
are separated by a deconfinement transition of vison con- 
densation, FL*pf-FLpf, corresponding to a proliferation of 
the (7r,a/2) fractional defects, that is expected to be in 
the inverted Ising universality class. 

The variety of phases and transitions between them are 
summarized in the phase diagram, FigjSjand a flow-chart, 
Fig|4] This rich fluctuations-driven phase behavior con- 
trasts sharply with a direct LO-N transition (described 



by C/(l)x C/(l) Landau theory Hmft = r{\A^ 



Ai(|A+i4 
theory. 



|A_|4) + A2|A+|^|A_|-^) found in mean-field 



G. Fermionic excitations 



The subtlety of the LO and its descendent states is in 
addition to strongly fluctuating Goldstone modes, they 
exhibit gapless fermionic excitations that can be strongly 
coupled to the superfluid phase (j) and smectic phonon u 
Goldstone modes. A complete comprehensive treatment 
remains an open problem. However, a number of approx- 
imate numerical and analytical treatments has lead to a 
consistent picture fSTl I134I[T77] . 

At large imbalance, near hc2 where A modulation is 
weak, the fermionic spectrum can be approximated by 
various q branches of the form Eq.([7| for the FF state. 
These lead to a novel paired superfluid with a Fermi sur- 
face of Bogoliubov quasi-particles. 

Near hd, a periodic form of A(r) is more appro- 
priately treated as an array of domain-walls of width 
f^ |134l[T57] . rather than a single harmonic. Diagonal- 
izing the BdG equations in their presence leads to a 
band of midgap states as in Id where this can be 
done exactly [S71 11471 1178j corresponding to imbalanced 
fermionic atoms residing on the nodes of the solitons 
in A akin to the polvacetylene |18(3] . Consistent with 
the large imbalance regime this leads to a Fermi sur- 
face pockets 179 illustrated in Fig|9] Despite consider- 
able progress a fully self-consistent quantum-mechanical 
treatment is still missing. 



H. Open questions 

While a broad range of phenomena associated with the 
Larkin-Ovchinnikov state has been explored, many inter- 
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FIG. 9: An illustration of Fermi pockets (full curve) of 
the gapless Bogoliubov quasi-particles characteristic of the 
Larkin-Ovchinnikov ground state. The periodic array of 
domain- walls in A_lo(2), the associated wavevector q, and 
the Fermi-surface of the underlying normal state (dashed cir- 
cle) are also indicated. 

esting and difficult questions remain open. Probably the 
most urgent of these is the long-standing question of the 
range of energetic stability of the crystalline supercon- 
ductor. If the state is indeed stable over a sufficiently 
broad range of detuning and imbalance to be experi- 
mentally accessible, is its lowest energy form indeed the 
striped coUinear LO type? While for large atomic clouds 
and shallow traps (such that LDA remains valid) only a 
small deformation of the LO state near the boundaries is 
expected, for tighter traps a more detailed treatment of 
the trap is necessary, and may lead to a distinct global 
form of the LO state, such as the "onion" and "radial" 
structures. To address such questions undoubtedly re- 
quires numerical solutions in experiment-specific geome- 
tries. 

Furthermore, the nature of the (2d and 3d) tran- 
sition into the LO state at the lower-critical Zeeman 
field /i^i . (135| . and the extent to which it resembles a 
commensurate-incommensurate transition (as in ld (134i 
I152| ) remains an open question. More broadly, while 
a number of LO descendent states have been proposed 
their detailed phenomenology, stability to quantum and 
thermal fluctuations, as well the nature of the associated 
phase transitions remains wide open. Similarly to the 
LO state, these phases are expected to exhibit gapless 
fermionic excitations coupled to their Goldstone modes. 
Understanding the effects of these fermionic modes on 
the phases and the associated transitions remain an ex- 
tremely interesting and challenging problem. 

To summarize, I reviewed a wide range of fluctuation 



phenomena in a LO state, expected to be realizable in 
an imbalanced resonant Fermi gas. Combining a micro- 
scopic analysis with robust model-independent symme- 
try arguments predicts that the LO state in an isotropic 
trap is a gapless superfluid smectic liquid crystal. Conse- 
quently, the state is extremely sensitive to thermal fluctu- 
ations that destroy its long-range positional order even 
in three dimensions, replacing it by a quasi-long range 
order, characterized by power-law correlations akin to 
a system tuned to a critical point or two-dimensional 
xy-model systems. This exotic state also exhibits vor- 
tex fractionalization, where the basic superfluid vortex is 
half the strength of a vortex in a regular paired conden- 
sate, and is accompanied by half-dislocations in the LO 
smectic (layered) structure. 

Analysis of the fluctuation-driven disordering of the 
LO smectic predicts a rich variety of descendant quantum 
liquid states, such as the superfluid (SFpf) and Fermi liq- 
uid (FLn) nematics and the fractionalized nonsuperfluid 
states {FL* ) , that generically intervene between the LO 
state and the conventional BCS superfluid (at low pop- 
ulation imbalance) and a conventional Fermi liquid (at 
high population imbalance). This phenomenology has a 
rich variety of experimental implications, many of which 
await detailed analysis. 



III. P-WAVE RESONANT BOSE GAS 
A. Background 

A p-wave resonant Bose gas is another system that 
was recently predicted to exhibit quantum liquid crys- 
tal order. Its study was in part motivated by the inter- 
esting phenomenology found in s-wave resonant bosonic 
systems (e.g., in ^^Rb ^26]), that was demonstrated to 
exhibit a magnetic field-driven quantum Ising transition 
between molecular and atomic superfluids [27H291 1183) , 
contrasting with a smooth BEC-BCS crossover in bal- 
anced fermionic isotopes. 

Although instabilities intrinsic to resonant bosons f5Bl 
I159j challenge realization of such bosonic molecular con- 
densates, many features of the phase diagram are ex- 
pected to survive away from the resonance and/or re- 
fiected in the nonequilibrium phenomenology (before the 
onset of the instability) of a resonant Bose gas. Further- 
more, recent extension of an s-wave resonant Bose gas to 
an optical lattice [1841 1185) demonstrated the stabiliza- 
tion through a quantum Zeno mechanism proposed by 
Rempe |186j . that dates back to Bethe's )187) analysis of 
the triplet linewidth in hydrogen. 

The predictions gTHlSl [HI ttSi] for the s-wave case 
have been supported by recent density matrix renor- 
malization group (188) . exact diagonalization [189 , and 
quantum Monte Carlo |190) studies. I am enthusiastic 
about similar progress in the substantially richer p-wave 
case, that I discuss below. 

Experimental realization of two-species degenerate 



15 



Bose gas with a p-wave Feshbach resonant inter-species 
interaction in ^^Rb - ^^Rb mixtures |159j provides a 
direct motivation for the theoretical studies discussed 
below[551 [55] , Related studies of Bose condensation in 
p-(and higher) bands in optical lattices have also been 
carried out[6Tl|62], but are only tangentially relevant to 
the present focus of quantum liquid crystals. 



B. Summary 

The phenomenology of a balanced two-component p- 
wave resonant Bose gas is summarized by a temperature- 
detuning phase diagram. As illustrated in Fig. [TO] at 




FIG. 10: Schematic temperature-detuning phase diagram for 
a balanced two-species p-wave resonant Bose gas. As il- 
lustrated, it exhibits atomic (ASF), molecular (MSF), and 
atomic-molecular (AMSF) superfluid phases. The novel 
AMSF state is characterized by a p-wave, molecular and a 
finite-momentum Q (see Fig. Ill atomic superfluidity. 



low temperature such two-component Bose gas generi- 
cally exhibits three classes of superfluid phases, atomic 
(ASF), molecular (MSF) and atomic-molecular (AMSF) 
condensates. 




FIG. 11: The momentum Qiiy) characteristic of the AMSF 
(polar) state, ranging between zero and the p-wave FR width- 
dependent value. 

The ASF appears at a large positive detuning (weak FR 
attraction) and low temperature, where one of the three 



combinations (ASFi, ASF2, ASF12) of the ^^Rh and 
*^Rb atoms are Bose-condensed into a conventional, uni- 
form superfluid, and the p-wave ^^Rb-®''Rb molecules are 
energetically costly and therefore appear only as gapped 
excitations. 

In the complementary regime of a large negative de- 
tuning, the attraction between two flavors of atoms is 
sufficiently strong so as to bind them into a tight p-wave 
hetero-molecules (e.g., ^^Rb-*^Rb molecule), which at 
low temperature condense into a p-wave molecular super- 
fluid isomorphic to a spinor-1 condensate [60, 191 200|. 
The latter is known to come in two forms, thereby pre- 
dicting the 4=0 "polar" (MSFp) and 4 = ±1 "ferro- 
magnetic" (MSFfm) molecular p-wave superfluid phases, 
with their relative stability determined by the ratio ag /a2 
of molecular spin-0 (aq) to molecular spin-2 (02) scatter- 
ing lengths. 

Besides these fairly conventional uniform atomic and 
molecular BEC's, for intermediate detuning around a 
unitary point the gas is predicted to exhibits novel 
AMSFp and AMSFfm phases, characterized by a non- 
zero momentum HQ atomic condensate j53[ I59( 161) , that 
is a superposition of two atomic species. The momentum 



Q 



< 



(40) 



is tunable via the FR detuning, primarily entering 
through the molecular condensate Umi^) ^ 1/^^ den- 
sity, and sensitive to the FR width 7p [3]. In addition 
to exhibiting an off-diagonal long-range order (ODLRO) 
of an ordinary superfluid the two AMSFp, states (dis- 
tinguished by the polar versus ferromagnetic nature of 
their p-wave molecular condensates) spontaneously par- 
tially break orientational and translational symmetries, 
akin to polar and smectic liquid crystals p8j and the pu- 
tative Fulde-Ferrell-Larkin-Ovchinnikov states of imbal- 
anced paired fcrmions P [5115111 ICT 1551 IMl 1711 [751151 fTUO] 
discussed in the flrst part of this review. This state is a 
finite momentum, HQ spinor superfluid, akin to (but dis- 
tinct from) a supersolid |125H128j . 

The physical picture behind such a finite-momentum 
AMSF formation is quite clear and is illustrated as a 
cartoon in AMSF phase of Fig. [TOj At intermediate de- 
tuning, where atomic gap closes within the MSF state, 
p-wave molecules decay via FR into a pair of atoms. 




FIG. 12: A cartoon of a p- wave molecule decaying into two op- 
positely moving two species of atoms, illustrating a resonant 
mechanism for a finite momentum Q atomic superfluidity (in- 
dicated by wavy lines) in the AMSF phase. 
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which (due to the p-wave nature of the molecules) are 
necessarily created at finite and opposite momenta, ±k, 
and therefore at low temperature form a finite mo- 
mentum atomic condensate, AMSF. The energetic cost 
(^ Q^/2to) of a finite momentum atomic condensation 
is balanced by the lowering of the energy (^ aQ^rim) 
through FR hybridization between closed-channel p-wave 
molecule and open-channel pair of atoms that is only pos- 



sible at finite atomic momentum, giving Q in Eq. ( 40 ) 



As illustrated in Fig. 13 in the polar AMSFp state, Q 



aligns along the quantization axis along which the molec- 
ular condensate has a zero projection of its internal I = \ 
angular momentum. For the case of the ferromagnetic 



£ = h X ifi 





FIG. 14: Schematic of the AMSFfm ferromagnetic state. The 
thick arrow indicates the atomic condensate momentum Q, 
lying in the plane transverse to the quantization axis £, along 
which the projection of the molecular internal orbital angular 
momentum is £z = +1- 



single-particle excitations are "squeezed" by the presence 
of the molecular condensate, offering a mechanism to re- 
alize atomic squeezed states 



FIG. 13: Schematic of the AMSFp polar state. The thick the MSF-AMSF transition at 



that can be measured 
by interference experiments, similar to those reported in 
Ref. 12021 The low-energy nature of these single-atom ex- 
citations is guaranteed by the vanishing of the gap at 

MSFp,f„-AMSFp,f„ 



with 



arrow indicates the atomic condensate momentum Q and the 
h arrow denotes the quantization axis along which the pro- 
jection of molecular internal orbital angular momentum van- 
ishes. 

AMSFfm state, Q lies in the otherwise isotropic plane, 
transverse to the p-wave molecular condensate axis, as 



illustrated in Fig. 14 



As any neutral superfluid, ASF, MSF, and AMSF are 
each characterized by Bogoliubov modes, with long wave- 
length acoustic "sound" dispersions 



<sf(^c) - 0. 

In addition to conventional Bogoliubov modes the 
AMSF exhibits a Goldstone mode corresponding to the 
fluctuations of a relative phase between the two atomic 
condensate components. A spatially periodic collinear 
AMSF state exhibits a smectic-like anisotropic phonon 
mode, akin to striped FFLO states [SH |66l 1^ SIH] dis- 
cussed in the first part of this review. 



(41) 



where (with a =ASFi,2,i2, MSFpjm, AMSFp^f^) are 
the associate d sound s peeds with standard Bogoliubov 
form Ca ~ \/ gcTia l^m. In each of these SF states one 
Bogoliubov mode (and only one in the ASF^ states) cor- 
responds to the overall condensate phase fiuctuations. 
In addition, the MSFp exhibits two degenerate "trans- 
verse" Bogoliubov orientational acoustic modes. The 
MSFfni is also additionally characterized by one "ferro- 
magnetic" spin-wave mode, i<;^^^f"> ^ fc2 ^nd one gapped 
mode, consistent with the characteristics of a conven- 
tional spinor-1 condensate |194[ 1195] . 

Because MSFp fm are paired molecular superfluids, 
they also exhibit gapped single atom-like quasiparticles 
(akin to Bogoliubov excitations in a fermionic paired BCS 
state), that do not carry a definite atom number. These 



WAMSFp(k) = ^{Bkl+Kk\)/x-, (42a) 

^AMSF..„(k) = ^{Bkl + k^{K,kl + Kykl))lx-, 

(42b) 

and an orientational mode wjl^^^, associated with orienta- 
tional symmetry breaking in AMSFfm 



I 

X+m 



k, 



Ijk^ [Bkl + k^{K^kl + Kykl, 



JX-k^ 



K^k^ 



(43a) 
(43b) 



where B, K^s, J, k, and x are compressional and bending 
moduli characterizing the phases. 
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C. Microscopic model of a balanced p-wave 
resonant Bose gas 

As required by bosonic statistics a p-wave reso- 
nance takes place between two distinguishable bosonic 

atoms created by i/'J(r) = (^^j\{r) , ipllr)^ , (e-g-, ^^Rb, 

^''Rb) |159| . interacting through a closed-channel £ — 
1 hetero-molecule created by a vector field operator 
(f>^{r) = {(j)l.,(j)^y,(f>l). The corresponding Hamiltonian is 
given by, 



bg 



(44) 



h.c. 



where = -jk^^ - A*(t, = -j^V^ - ^i,n, with 
molecular chemical potential /im = /^i + M2 — ad- 
justable by a magnetic field dependent detuning v, latter 
the rest energy of the closed-channel molecule relative 
to a pair of open-channel atoms. A generalization to 
number and mass imbalanced mixtures (as studied for 



fermionic atoms [1 O [SHI [Ml IHl IZi 111 [100] remain to 
be explored [204] , 

For simplicity, above form is specialized to a rotation- 
ally invariant FR interaction, with w and a independent 
of the molecular component i. This is an approximation 
for the system of interest, ^^Rb-^^Rb mixture, where in- 
deed the p-wave FR around B = 257.8 Gauss [159) is 
split into a doublet by approximately AS = 0.6 Gauss, 
similar to the fermionic case of ''^K [31 [H [H [203] ■ A 
more realistic, richer case remains to be explored. 

The background (non-resonant) interaction density 
Hbg is the short-scale two-body interaction between spin 
and number densities characterized by atomic scattering 
lengths ai, ai2 as well as molecular scattering lengths. 
A miscibility of a two-component atomic gas requires a 
condition on the corresponding atomic s-wave scattering 
lengths aia2 > a^2 ^205' , which may be problematic for 
the case of *^Rb-*^Rb due to the negative background 
scattering length of *^Rb. 

The corresponding imaginary time (r) coherent state 
Lagrangian density is given by 



V2 

2m 

+ Ai2iV'inv'2p 



y2 

4m 

(|^lP + |V'2p)|</'P 



91 



2 , 92 



+ -^i4>* ■ [V'l(-jV)V'2 - V'2(~«V)V'i] + C.C.) , 



(45) 



where A^s and ^^s are the atomic and molecular s-wave 
pseudopotentials. Closely related models also arise in 
completely distinct physically contexts. These include 
quantum magnets that exhibit incommensurate spin liq- 
uids states |176j and bosonic atoms in the presence of 
spin-orbit interactions [54] . 

The low-energy two-atom vacuum scattering in the 
above two-channel model can be computed exactly. It 
faithfully captures all the features of the low-energy p- 
wave resonant and s-wave nonresonant scattering phe- 
nomenology of the ^^Rb-^^Rb p-wave Feshbach-resonant 
mixture |159| . By matching experiments in the dilute 
limit one can fix the model's key parameters [51 155] . The 
analysis at nonzero, balanced atomic densities leads to 
the predictions summarized above. All predictions are 
qualitatively robust, and can furthermore be made quan- 
titatively accurate in a narrow-resonance limit. 



of the Lagrangian, ( 45 1 , supplemented by symmetry ar- 



guments and the analysis of fluctuations about the or- 
dered state. 

In the absence of periodic or disorder potential at suf- 
ficiently low temperatures a bosonic gas is always a su- 
perfluid, that in three dimensions exhibits Bose-Einstein 
condensation, characterized by scalar atomic, ^'o- and/or 
3- vector molecular, # complex order parameters. Thus, 
at low T the gas exhibits three classes of SF phases: 

1. Atomic Supcrfiuid (ASF), 7^ and * = 

2. Molecular Superfluid (MSF), = and * 7^ 

3. Atomic Molecular Superfluid (AMSF), 7^ and 
* 7^ 

with the transition between them driven by the magnetic 
field-dependent detuning, ly. 



D. Phases, symmetries and Goldstone modes 



1. Atomic superflmds: ASFi, ASF2, ASF12 



Qualitative features of the phase diagram for this sys- 
tem can be mapped out through a mean-field treatment 



For large positive detuning i>, closed-channel molecules 
are gapped and the ground state is a molecular vac- 



18 



uum, and a zero-momentum atomic condensate ASF. 
The latter itself comes in three forms: (i) ASFi with 
*i 7^ 0,^-2 = 0, (ii) ASF2 with *i = 0,*2 7^ 0, (iii) 
ASF12 with ^1 ^ 0, ^2 7^ 0, separated by continuous 
phase transitions. For a balanced mixture /2i — ji2, the 
system exhibits a direct N-ASF12 transition through a 
tetracritical point, fli = — 0, that is believed to be 
in the decoupled universality class [206 - 208 . All other 
transitions (N-ASFi, N-ASF2, and ASF^-ASFiz) are in 
the XY universality class, breaking associated U{1) sym- 
metries. The phase boundaries and the values of the 
atomic condensate order parameters can be straightfor- 
wardly computed within mean- field theory [S31 [52] , but 



are modified by fluctuations |207| r208j . 

Within ASF phases, the spectrum of fluctuations can 
be straightforwardly computed by a Bogoliubov diago- 
nalization of coupled atomic and molecular excitations, 
with details depending on which of the three possible 
ASF phases is studied. In general these states exhibit one 
Bogoliubov sound mode per broken atomic U(l) symme- 
try, with one Goldstone mode in ASFi and ASF2 phases 
and two in ASFi2i59]. 



2. Molecular super-fluids: MSFp and MSFfm 

In the opposite limit of a large negative detuning, 
atoms are gapped, tightly bound into heteromolecules, 
that at low temperature condense into a p-wave molecu- 
lar superfluid, MSF. In this regime of atomic vacuum, the 
gas reduces to that of interacting p-wave molecules, iso- 
morphic to that of the extensively studied F — 1 spinor 
condensate [5D1 119111200] , with the hyperfine spin F here 
replaced by the orbital £ — I angular momentum of two 
constituent atoms. 

Like F ^ I spinor condensates, the p-wave molecular 
superfluid, MSF exhibits two distinct phases depending 
on t he s ign of the renormalized interaction coupling 32 in 



(m) 



Eq. ( 45 1 , or equivalently the sign of the difference Oq 

aj™"* of the molecular L = and L = 2 channels s-wave 
scattering lengths. 

Polar molecular superfluid, MSFp: 

For 52 < the ground state is the so-called "po- 
lar" molecular superfluid, MSFp, characterized by a 
(collinear) order parameter # = ^pe^'^fi, with h a real 
unit vector, tp a (real) phase, and $p a (real) order- 
parameter amplitude, with the state corresponding to 

= projection of the internal molecular orbital an- 
gular momentum along n. MSFp clearly spontaneously 
breaks rotational symmetry by its choice of the £z = 
quantization axis n, and the global gauge symmetry, cor- 
responding to a total atom number conservation. The 
low-energy order parameter manifold that characterizes 
MSFp is given by the coset space (C/(l) (g) S2)/^2, ad- 
mitting half-integer "charge" vortices |200j akin to (but 
distinct from) the s-wave MSF 



wav|59|. One corresponds to breaking of the global atom 
number conservation and two associated with breaking 
of rotational 0(3) symmetry [5^ I194j . 

Ferromagnetic molecular superfluid, MSFfm- 
Alternatively, for 172 > the ground state is the "fer- 
romagnetic" molecular superfluid, MSFfm, characterized 
by an (coplanar) order parameter $ = ^^{n±irh), with 

n, rh, £ = n X rh a real orthonormal triad, <I>tm a real 
amplitude, and the state corresponds to £z = ±1 projec- 
tion of the internal molecular orbital angular momentum 
along the £ axis. MSFfm spontaneously breaks the time 
reversal, the 0(3) rotational and the global gauge Un{^) 
symmetries, latter corresponding to a total atom number 
N conservation. Inside MSFfm the low-energy order pa- 
rameter manifold is that of the 0(3) — SU{2)/Z2 group, 
corresponding to orientations of the orthonormal triad 
n, rh, £. 

As its hyperfine spinor-condensate cousin, the ferro- 
magnetic MSFfm exhibits (only) two gapless Goldstone 
modes, one linear (oc k) conventional Bogoliubov mode 
associates with the broken global U{1) gauge symmetry, 
and another quadratic (cx fc^) corresponding to the ferro- 
magnetic order, with associated spin-waves |194) reflect- 
ing the precessional FM dynamics. This is despite the 
three-dimensionality of its 50(3) coset space and can be 
traced back to the fact that the two components of the 
spinor are canonically conjugate and as a result combine 
into a single low-frequency mode. 



MSFta 




pMSFta 



^MSFr, 



FIG. 15: Excitation spectrum for the ferromagnetic molecu- 
lar superfluid, MSFfm- The doubly-degenerate atomic spec- 
trum (thin curves) exhibits a minimum gap at nonzero k, a 
precursor of finite momentum atomic condensation into the 
AMSFfm. The molecular spectrum (thick curves), consists 
of a longitudinal gapless quadratic ferromagnetic spin-wave 
mode (lowest), a Bogoliubov sound mode and a quadratic 
gapped mode. 



As illustrated in Fig. (15), while atomic excitations are 



gapped (confined into molecules), because of the ^i-wave 
FR coupling in ( 45 1 the minimum in the atomic spectrum 



The polar MSFp state exhibits three gapless 
Bogoliubov-like modes that are calculated in a standard 



occurs at a nonzero momentum, Qp, a precursor of the 
Bose-condensation transition into a nonzero momentum 
AMSF state, when this atomic gap closes. 
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3. Atomic-molecular superfluids: AMSFp and AMSFfm 

As detuning is further increased from a large negative 
value of the MSFp^tm phases, for intermediate v the gap 
to atomic excitations decreases, closing at a critical value 
of j/MSF-AMSF g^j-^j leading to atomic Bose-condensation 
into the corresponding AMSFp and AMSF/m states. It 
is clear from the linear momentum dependence of the p- 
wave FR coupling in the Hamiltonian (44), that, quite 



generally, the MSF-AMSF transition is robustly into a 
nonzero-momentum atomic condensate, with k — Q, set 
by a balance of the p-wave FR hybridization and the 
atomic kinetic energies. 

As with other (partially) crystalline states of mat- 
ter [49H5T] . the nature of the resulting AMSF states de- 
pends on the symmetry of the crystalline order, encoded 
into the atomic condensate order parameter 



*2(r) 




(46) 



via the associated set of the reciprocal lattice vectors, 
at which the condensation takes place. Determined by a 
detailed nature of interactions and fluctuations, typically 
the nature of crystalline order is challenging to deduce. 

The nature and symmetries of these AMSF states fur- 
thermore qualitatively depends on the parent MSF, with 
the polar AMSFp and the ferromagnetic AMSFf,„ as two 
possibilities. In addition to the symmetries already bro- 
ken in its MSF parent, by virtue of atomic condensa- 
tion the AMSF state breaks the remaining C/AAr(l) global 
gauge symmetry associated with the conservation of the 
difference in atom species number, AiV. Other symme- 
tries that it breaks depend on the detailed structure of 
the AMSFfm^p states. 

Polar atomic-molecular superfluid, AMSFp 



atomic-molecular superfluid AMSFp 
the polar molecular state, MSFp. As 



The polar 
emerges from 

illustrated in Fig 13 and is clear from the form of (45) 
in the AMSFp the finite momentum atomic condensate 
orders with Q locked along the molecular condensate 
field # = <I>pe*'^n. Hence, quite generically it is a 
Fulde-Ferrell like ^50] collinear plane-wave state, that 
spontaneously breaks the time reversal symmetry, but 
not any additional spatial symmetries. A schematic 
phase diagram for this polar class of states is illustrated 
in FigjlG] The phase boundaries corresponding to the 



MSFp 
given by 



AMSFp and the AMSFp 



ASF transitions are 



,MSF„-AMSF„ 



,AMSFp-ASF 



= - (ffi + 52 ~ 2gQ„i) n„, (47a) 

~ - ^ (.91 +92- 2gam) n, (47b) 

= {2X~ ga,n)na, (47c) 

w(2A-.9„„)n, (47d) 



The evolution of the order parameters across the phase 




$ = 0,^_ 7^ 

ASF 




FIG. 16: Mean field phase diagrams for polar phase as a func- 
tion of atomic and molecular chemical potentials, Ha, t^m, re- 
spectively. Ferromagnetic phase is similar but with different 
parameters, (a) For \{gi -\- (72) — gam > 0, all three super- 
fluid phases, ASF, AMSF, and MSF appear and are separated 
by continuous phase transitions (thick black lines), (b) For 
Hdi + 92) — gam < 0, AMSF is unstable, and the ASF and 
MSF are separated by a first-order transition (hatched double 
line) . 



Despite a nonzero momentum (Q) of its condensate, 
in its minimal form such AMSFp state is a nematic (not 
smectic) superfluid that is anisotropic but spatially ho- 
mogeneous. A detailed analysis [SS] shows that Goldstone 
modes in this phase are governed by a Lagrangian density 



5L^ 



2 m 

I TO 



(48) 



K 



(Vi^?_) 



diagram is illustrated in Fig 17 



where the compressibilities x± ^'^d other couplings are 
functions of the microscopic parameters that have been 
computed in weakly interacting limit [59]. Thus, despite 
a translationally invariant nature of the AMSFp phase, 
in addition to the conventional (linear in fc) Bogoluibov 
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MSF, 



Vc2 ASF 



FIG. 17: Atomic (thick) and molecular (thin) order param- 
eters versus the FR detuning v for the polar phase, with 



MSF„-AMSF, 



and Vci = v\ 



AMSFn-ASF 



mode, 6*+, akin to the FF state[Sni[Mlin3[nH| it exhibits 
a smeciic-like Goldstone mode, with dispersion given 
in Eq.(|42a|[5^. 



Ferromagnetic atomic-molecular superfluid, AMSFfm 
Similarly, upon increase of the detuning, a finite- 



momentum atomic condensation from the MSFfm 
leads to the ferromagnetic atomic-molecular superfluid, 
AMSFfni- In this state, a p- wave Feshbach resonant inter- 
action leads to the energetic preference for a transverse 
orientation of the atomic condensate momentum Q to 
the molecular quantization axis, £ = n x rh. Conse- 
quently, as illustrated in FigjMj the AMSFfm state spon- 
taneously breaks additional rotational symmetry of the 
uniaxial molecular state in the plane perpendicular to 
the molecular quantization axis £. In fact the nature of 
spatial order inside the n — rh plane is not necessarily 
of co-linear striped order, and generically admits a 2d 
crystalline condensate, nontrivially determined by inter- 
actions and fluctuations. Motivated by the polar version 
of AMSF, so far only a striped order has been inves- 
tigated, with general form remaining an open problem. 
For such collinear order, a mean-field analysis predicts 
a FF (rather than LO jM]) type AMSF/™ state in the 
h — rh plane. Thus, the AMSFfm state is a biaxial ne- 
matic superfluid, defined by Q and £ axes. 

The corresponding Goldstone modes Lagrangian den- 
sity is given by: 



+ iKdye.dr-i + \j{v-if, 



^f + -Ky{vdye,f 



(49) 



where the Goldstone mode 7 describes one additional 
fluctuating angle of the rh axis outside of the h-rh plane, 
into the £ axis. 

The biaxiality is expected and arises due to the vector 
p-wave order, characterized by a spinor $/m, with the 
quantization axis, £. The finite angular momentum, — 
±1 along £ distinguishes AMSFfm from AMSFp and leads 
to this additional Goldstone mode 7. 

A straightforward diagonalization of the above La- 
grangian leads to dispersions for three Goldstone modes 
inside the AMSFfm state: 



^[Bkl + k^{K^kl+Kykl)]/x^ 
ljk-^[Bkl + k-^{K,kl + Kykl)] 



JX-k^ 



(50a) 
(50b) 

(50c) 



The anisotropic Ci;J^(k) dispersion corresponds to the 
ferromagnetic spin-waves in the plane of atomic conden- 
sate phase-fronts ( "smectic layers" ) of the p-wave atomic- 
molecular condensate, AMSFfm, reducing to the disper- 
sion of MSFfm for a vanishing smectic order, with B = 0. 



E. Phase transitions 

Phase transitions from the high temperature normal 
state into the ASF and MSF states have been well- 
explored in the context of conventional and spinor su- 
perfluids. 

The quantum MSF - AMSF phase transitions are more 
interesting and require a beyond mean-field treatment. 
In large part the complexity is associated with the non- 
trivial coupling of the Goldstone modes of the MSF phase 
to the finite-momentum order parameter condensing in 
the AMSF state. 

Starting with the coherent-state microscopic La- 
grangian, specializing to the MSFp ordered state, and 
integrating out massive modes, one obtains the effective 
Lagrangian density 

£p ^e^^\drij-f + ^1 (-iV - amyfp:,^5h) tp^f 



2 2g2 



2 , ^™0(^.^2 



4m 
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(V¥')^ 



(51) 



2g 4m 
that describes the quantum MSFp-AMSFp transition 
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The critical field ■0- 



^("01,(5 + "02 -q) is the finite- 
momentum atomic condensate order parameter, n the 
^2=0 quantization axis and (f) the superfluid phase, 
defining the MSFp state. 

Thus, as can also be argued on symmetry grounds, 
the zero-temperature transition is indeed described by a 
quantum ((d -f- l)-dimensional) de Gennes model (also 
known as the Ginzburg-Landau or Abelian-Higg's mod- 
els) |48j , akin to that for a normal-to-superconductor and 
nematic-to-smectic-A transitions. The role of the ne- 
matic director (gauge-field) is played by the = quan- 
tization axis of the p-wave molecular condensate. Based 
on the extensive work for these systems [2091 I210j . in 
three (spatial) dimensions (d = 3) the effective gauge- 
field fluctuations drive this transition first-order. 

Similarly, the quantum MSFfm-AMSFtm transition is 
described by the effective low-energy Lagrangian density 



C 



fm 



1 
2m 



-iV 

PmQ 

8m 



-on V- 



V2 



8fn 



(52) 



It is distinguished from the polar case by the additional 
biaxial order whose fluctuations are characterized by rh. 



F. Topological defects 



in the ASF^ states and well explored defects in the spin-1 
AMSFp condensates, the AMSFp exhibit interest- 
ing combination of vortex, dislocation and domain-wall 
defects similar to those discussed in the context of FFLO 
and other striped superfluids. While some preliminary 
studies have been carried out [59 , their detailed analysis, 
response to rotation, and realization remain subjects of 
future research. 



IV. SUMMARY AND CONCLUSION 

In this review I discussed two recently studied quan- 
tum liquid crystal orders, proposed to be realizable in 
s-wave resonant fermionic and p-wave resonant bosonic 
atomic gases. The key ingredient is the interaction- 
driven superfluid condensation at a nonzero momentum, 
that in an isotropic trap spontaneously partially breaks 
orientational and in some cases translational symmetry. 
The main state that generically appears at low temper- 
ature is a quantum superfluid smectic. A rich variety of 
other phases results from partially disordering this state 
by unbinding various combinations of topological defects 
associated with compact nature of the broken symme- 
tries. An interesting interplay ensues between Goldstone 
modes, associated spontaneously broken gauge and spa- 
tial symmetries and the fermionic and bosonic single- 
particle excitations. Although a general formulation of 
such states is credibly established many interesting open 
questions remain. 



In addition the smooth low-energy configurations of 
these Goldstone modes, the compact nature of the or- 
der parameters characterizing the ASF, MSF and AMSF 
states admits a rich variety of topological defects. These 
singular excitations are crucial to a complete character- 
ization of the states and their disordering, particularly 
in the case of non-meanfield (e.g., partially disordered) 
states that are not uniquely characterized by a Landau 
order parameter. In addition to the conventional vortices 
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